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Abstract 

The cross section for high-i^^ dijet production in photoproduction has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
81.8 pb~^. The events were required to have a virtuahty of the incoming pho- 
ton, Q^, of less than 1 GeV^ and a photon-proton center-of-mass energy in the 
range 142 < W^p < 293 GeV. Events were selected if at least two jets satis- 
fied the transverse-energy requirements of E}^^^ > 20 GeV and E}^^"^ > 15 GeV 
and pseudorapidity (with respect to the proton beam direction) requirements of 
— 1 < r/j'^*^'^ < 3, with at least one of the jets satisfying —1 < rf'^^ < 2.5. The 
measurements show sensitivity to the parton distributions in the photon and 
proton and to effects beyond next-to- leading order in QCD. Hence these data 
can be used to constrain further the parton densities in the proton and photon. 
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1 Introduction 



In photoproduction at HERA, a quasi-real photon emitted from the incoming positron^ 
colhdes with a parton from the incoming proton. The photoproduction of jets can be 
classified into two types of processes in leading-order (LO) Quantum Chromodynamics 
(QCD). In direct processes, the photon participates in the hard scatter via either boson- 
gluon fusion (see Fig. 1(a)) or QCD Compton scattering. The second class, resolved 
processes (see Fig. 1(b)), involves the photon acting as a source of quarks and gluons, 
with only a fraction of its momentum, x^, participating in the hard scatter. Measurements 
of jet cross sections in photoproduction [1-6] are sensitive to the structure of both the 
proton and the photon and thus provide input to global fits to determine their parton 
densities. 

There are three objectives of the measurement reported in this paper. Firstly, the analysis 
was designed to provide constraints on the parton density functions (PDFs) of the photon. 
Over the last two years there has been active research in the area of fitting photon PDFs 
and a number of new parameterizations have become available [7-9]. In two of these [7,8], 
fits were performed exclusively to photon structure function, F^, data; the other [9] 
also considered data from a previous dijet photoproduction analysis published by the 
ZEUS collaboration [4]. It is the purpose of this analysis to test the effectiveness of each 
parameterization at describing HERA photoproduction data. To this end, the present 
analysis was conducted at higher transverse energy relative to previous publications. It 
is expected that at these high transverse energies the predictions of next-to-leading-order 
(NLO) QCD calculations should describe the data well, have smaller uncertainties, and 
allow a more precise discrimination between the different parameterizations of the photon 
PDFs. The reduction in statistics associated with moving to higher transverse energies 
was in part compensated by the factor of two increase in luminosity, for this independent 
data sample, and the extension to higher pseudorapidity^ of the jet compared to the 
previous analysis [4]. 

Secondly, the present analysis was designed to provide constraints on the proton PDFs. 
Global fits to determine the proton PDFs continue to be a very active and important 
area of research. A common feature of these global fits is a large uncertainty in the gluon 
PDF for high values of Xp, the fractional momentum at which partons inside the proton 

^ In the following, the term "positron" denotes gcncrically both the electron (e~) and positron (e+). 
Unless explicitly stated, positron will be the term used to describe both particles. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity 
is defined as 7] = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 
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are probed. At such high values {xp > 0.1), the gluon PDF is poorly constrained and so 
attempts were made for the present investigation to measure cross sections which show 
particular sensitivity to these uncertainties. Recently, the ZEUS collaboration included 
jet data into fits for the proton PDFs [10]. 

Finally, the difference in azimuthal angle of two jets was considered, as in previous mea- 
surements of charm and prompt photon photoproduction [11, 12]. In LO QCD, the cross 
section as a function of the azimuthal difference would simply be a delta function located 
at TT radians. However, the presence of higher-order effects leads to extra jets in the final 
state and in values less than vr radians. The cross section is therefore directly sensitive 
to higher-order topologies and provides a test of NLO QCD and of Monte Carlo (MC) 
models with different implementations of parton-cascade algorithms. The data for charm 
photoproduction [11] demonstrated the inadequacy of NLO QCD, particularly when the 
azimuthal angle difference was significantly different from vr and for a sample of events 
enriched in resolved-photon processes. To investigate this inadequacy in a more inclusive 
way and with higher precision, such distributions were also measured. 



2 Definition of the cross section and variables 

Within the framework of perturbative QCD, the dijet positron-proton cross section, dagp, 
can be written as a convolution of the proton PDFs, fp, and photon PDFs, with the 
partonic hard cross section, d&ab, as 



where y = E^/ Ef, is the longitudinal momentum fraction of the almost-real photon emitted 
by the positron and the function /^/e is the flux of photons from the positron. The equation 
is a sum over all possible partons, a and h. In the case of the direct cross section, the 
photon PDF is replaced by a delta function at = 1. The scales of the process are the 
renormalization, yU/j, and factorization scales, ^p- 

To probe the structure of the photon, it is desirable to measure cross sections as functions 
of variables that are sensitive to the incoming parton momentum spectrum, such as the 
momentum fraction, x-^, at which partons inside the photon are probed. Since x^ is not 
directly measurable, it is necessary to define [1] an observable, x'^^, which is the fraction 
of the photon momentum participating in the production of the two highest transverse- 
energy jets (and is equal to x^ for partons in LO QCD), as: 




ab 
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where is the incident positron energy, E^^^^^ and E^^^'^ are the transverse energies and 
^jeti g^^^ j^jet2 pseudorapidities of the two jets in the laboratory frame {E}^^^ > E^r^^"^). 
At LO (see Fig. 1), direct processes have x'^^ = 1, while resolved processes have x"^^ < 1. 

For the proton, the observable x"^^ is similarly defined [1] as 



^obs ^ -^T *^ -t- -C/T^ e 



4" = .i^;^^ (3) 

where Ep is the incident proton energy. This observable is the fraction of the proton 
momentum participating in the production of the two highest-energy jets (and is equal 
to Xp for partons in LO QCD). 

Cross sections are presented as functions of x*^^, x*^^, Et, E}^'\ f] and |A</>jj|. The mean 
transverse energy of the two jets, Et, is given by 

Et = ^ \^ ■ (4) 
Similarly, the mean pseudorapidity of the two jets, f/, is given by 

^jetl _j_ „jet2 

n = — ^ — • (5) 

The absolute difference in azimuthal angle of the two jets, 0^°*^ and is given by 



|A0jj| = -(/^''^^i^ (6) 

The kinematic region for this study is defined as < 1 GeV^, where = 2E(.E'^{1 + 
cos^e) and E'^ and 6e are the energy and angle, respectively, of the scattered positron. 
The photon-proton center-of-mass energy, W^p = ^jAyEeEp, is required to be in the 
range 142 GeV to 293 GeV. Each event is required to have at least two jets reconstructed 
with the kx cluster algorithm [13] in its longitudinally invariant inclusive mode [14], with 
at least one jet having transverse energy greater than 20 GeV and another greater than 
15 GeV. The jets are required to satisfy —1 < ■/^J'^*^'^ < 3 with at least one jet lying in the 
range between —1 and 2.5. The upper bound of 3 units represents an extension of the 
pseudorapidity range by 0.6 units in the forward direction over the previous analysis [4], 
thereby increasing the sensitivity of the measurement to low-x-y and high-Xp processes. 
The cross sections for all distributions have been determined for regions enriched in direct- 
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and resolved-photon processes by requiring x'^^ to be greater than 0.75 or less than 0.75, 
respectively. 

One of the goals of the present investigation is to provide data that constrain the gluon 
PDF in the proton, which exhibits large uncertainties at values of Xp > 0.1. A study 
was performed [15] by considering the x"^^ cross section in different kinematic regions, 
varying the cuts on the jet transverse energies and pseudorapidities as well as on x"^^. 
This allowed the determination of kinematic regions in which the cross section was large 
enough to be measured and in which the uncertainties on the cross section that arise 
due to those of the gluon PDF were largest. These cross sections will be referred to as 
"optimized" cross sections and are those which have the largest uncertainty from the gluon 
PDF; in total eight cross sections were measured (four direct enriched and four resolved 
enriched). The PDF sets chosen to conduct the optimization study were the ZEUS-S [16] 
and ZEUS- JETS [10] PDF sets. The kinematic regions of the cross sections are defined 
in Table 1, where the W^p and requirements are as above. 



3 Experimental conditions 

The data were collected during the 1998-2000 running periods, where HERA operated 
with protons of energy Ep = 920 GeV and electrons or positrons of energy = 27.5 
GeV. During 1998 and the first half of 1999, a sample of electron data corresponding 
to an integrated luminosity of 16.7 ± 0.3 pb~^ was collected. The remaining data up to 
the year 2000 were taken using positrons and correspond to an integrated luminosity of 
65.1 ± 1.5 pb~^. The results presented here are therefore based on a total integrated 
luminosity of 81.8 ± 1.8 pb~^. A detailed description of the ZEUS detector can be found 
elsewhere [17, 18]. A brief outline of the components that are most relevant for this 
analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [19], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering the 
polar-angle region 15° < 9 < 164°. The transverse- momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pr, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [20] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections 
(HAG). The smallest subdivision of the calorimeter is called a cell. The CAL energy res- 
olutions, as measured under test-beam conditions, are a{E)/E = 0.18/ ^/E for electrons 
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and a{E)/E = 0.35//E for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > ejp, 
where the photon was measured in a lead-scintillator calorimeter [21] placed in the HERA 
tunnel at Z = —107 m. 

4 Monte Carlo models 

The acceptance and the effects of detector response were determined using samples of 
simulated events. The programs Herwig 6.505 [22] and Pythia 6.221 [23], which im- 
plement the leading-order matrix elements, followed by parton showers and hadronization, 
were used. The Herwig and Pythia generators differ in the details of the implemen- 
tation of the leading-logarithmic parton-shower models and hence are also compared to 
the measured cross-section da/d\A(j^^\. The MC programs also use different hadronization 
models: Herwig uses the cluster model [24] and Pythia uses the Lund string model [25]. 
Direct and resolved events were generated separately. For the purposes of correction, the 
relative contribution of direct and resolved events was fitted to the data. For all gener- 
ated events, the ZEUS detector response was simulated in detail using a program based 
on Geant 3.13 [26]. 

For both MC programs, the CTEQ5L [27] and GRV-LO [28] proton and photon PDFs, 
respectively, were used. The p™™ for the outgoing partons from the hard scatter was set to 

4 GeV. For the generation of resolved photon events, the default multiparton interaction 
models [29,30] were used. A comparably reasonable description of the raw data kinematic 
distributions was observed with both Herwig and Pythia MC simulations. 

5 NLO QCD calculations 

The calculation for jet photoproduction used is that of Frixione and Ridolfi [31,32], which 
employs the subtraction method [33] for dealing with the coUinear and infra-red divergen- 
cies. The number of fiavors was set to 5 and the renormalization and factorization scales 
were both set to {E^^^^"^), which is half the sum of the transverse energies of the final-state 
partons. The parton densities in the proton were parameterized using CTEQ5M1 [27]; 
the value as{Mz) = 0.118 used therein was adopted for the central prediction. 

The following parameterizations of the photon PDFs were used: Cornet et al. (CJK) [7], 
Aurenche et al. (AFG04) [8], Slominski et al. (SAL) [9], Gliick et al. (GRV-HO) [28] and 
a previous set of PDFs from Aurenche et al. (AFG) [34]. The three new PDFs [7-9] use 
all available data on F2 from the LEP experiments. The data are of higher precision and 
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cover a wider region of phase space, reaching lower in and higher in the momentum of 
the exchanged photon, compared to the data used in the AFG and GRV-HO parameteri- 
zations. The parameterization from CJK uses a more careful treatment of heavy quarks, 
whereas that from SAL also considers previous dijet photoproduction data from ZEUS [4]. 
The most striking difference between the resulting PDFs is that CJK has a more rapid 
rise of the gluon density at low x^. 

The NLO QCD predictions were corrected for hadronization effects using a bin-by-bin 
procedure according to da = da^^^ ■ Chad, where da^^'^ is the cross section for partons 
in the final state of the NLO calculation. The hadronization correction factor, Chad, was 
defined as the ratio of the dijet cross sections after and before the hadronization process, 
Chad = (i(T^^''°'^Y'^'^MC*°°^- '^^^ value of Chad was taken as the mean of the ratios obtained 
using the Herwig and Pythia predictions. The hadronization correction was generally 
below 10% in each bin. 

Several sources of theoretical uncertainty were investigated, which are given below with 
their typical size, 

• the renormalization scale was changed to 2^"'^ [10]. This led to an uncertainty 
of t(10 - 20)%; 

• the factorization scale was changed to 2^^'^ ■ {E^^'^^"^) [10]. This led to an uncertainty 
of ±(5 - 10)%; 

• the value of as was changed by ±0.001, the uncertainty on the world average [35], 
by using the CTEQ4 PDFs for ^.(M^) = 0.113, 0.116 and 0.119 and interpolating 
accordingly. This led to an uncertainty of about ±2%; 

• the uncertainty in the hadronization correction was estimated as half the spread be- 
tween the two MC correction factors. This led to an uncertainty of generally less than 
±5%. 

The above four uncertainties were added in quadrature and are displayed on the figures 
as the shaded band around the central prediction. The size of these uncertainties is also 
shown as a function of Et, x"^^ and x*^^ in Fig. 2. The uncertainty from changing the 
renormalization scale is dominant. It should be noted that here the renormalization and 
factorization scales were varied independently by factors of 2^°'^ and the resulting changes 
were added in quadrature as in the determination of the ZEUS- JETS PDF [10]. The result 
of this procedure leads to an uncertainty which is approximately the same as varying both 
simultaneously by 2^^ as has been done previously [4]. 

Other uncertainties which were considered are: 

• the uncertainties in determining the proton PDFs were assessed by using the ZEUS- 
JETS PDF uncertainties propagated from the experimental uncertainties of the fitted 
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data. This led to an uncertainty of ±(5 — 10)%; 

• the uncertainties in determining the photon PDFs were assessed by using sets from 
different authors. Differences of generally less than 25% were observed between the 
AFG, AFG04, SAL and GRV sets. However, the predictions based on CJK were up 
to 70% higher than those based on the other four. 

These uncertainties were not added in quadrature with the others, but examples of their 
size are given in Fig. 2. Differences between the two photon PDFs, CJK and AFG04, are 
concentrated at low x'^^ and low Ex', the low x"^^ region is most sensitive to the gluon 
distribution in the photon, which increases more rapidly for CJK as shown in Fig. 3. At 
lowest 0;°'^'', the fraction of the cross section arising from the gluon distribution in the 
photon is 66% for CJK. The uncertainty on the proton PDF increases with increasing 
Ex and x'^^ and is sometimes, as seen in Fig. 2(c), as large as the other combined un- 
certainties. The fraction of the cross section arising from the gluon distribution in the 
proton is about 50% for the lower Ex and xf'^ values considered, but decreases to below 
20% for high values. However, the uncertainty on the gluon dominates the proton PDF 
uncertainty in most of the kinematic region investigated. 



6 Event selection 

A three-level trigger system was used to select events online [2,18,36]. At the third level, a 
cone algorithm was applied to the CAL cells and jets were reconstructed using the energies 
and positions of these cells. Events with at least one jet, which satisfied the requirements 
that the transverse energy exceeded 10 GeV and the pseudorapidity was less than 2.5, 
were accepted. Dijet events in photoproduction were then selected offline by using the 
following procedures and cuts designed to remove sources of background: 

• to remove background due to proton beam-gas interactions and cosmic-ray showers, 
the longitudinal position of the reconstructed vertex was required to be in the range 
I ^vertex I < 40 cm; 

• a cut on the ratio of the number of tracks assigned to the primary vertex to the 
total number of tracks, N^^^/N^rk > 0.1, was also imposed to remove beam- related 
background, which have values of this ratio typically below 0.1; 

• to remove background due to charged current deep inelastic scattering (DIS) and 
cosmic-ray showers, events were required to have a relative transverse momentum of 
px/y/Ex < 1.5 VCeV, where px and Ex are, respectively, the measured transverse 
momentum and transverse energy of the event; 



7 



• neutral current (NC) DIS events with a scattered positron candidate in the CAL were 
removed by cutting [1] on the inelasticity, y, which is estimated from the energy, E'^, 
and polar angle, 6'^, of the scattered positron candidate using j/e = 1 — — cos6'g). 
Events were rejected if i/e < 0.7; 

• the requirement 0.15 < yjB < 0.7 was imposed, where yjB is the estimator of y 
measured from the CAL energy deposits according to the Jacquet-Blondel method [37]. 
The upper cut removed NC DIS events where the positron was not identified and which 
therefore have a value of yjB close to 1. The lower cut removed proton beam-gas events 
which typically have a low value of yjB] 

• the /cr-clustering algorithm was applied to the CAL energy deposits. The transverse 
energies of the jets were corrected [3,4,38] in order to compensate for energy losses in 
inactive material in front of the CAL. Events were selected in which at least two jets 
were found with E^^" > 20GeV, E^^^^ > 15GeV and -1 < V^"'^ < 3, with at least 
one jet lying in the range between —1 and 2.5. In this region, the resolution of the jet 
transverse energy was about 10%. 



7 Data correction and systematics 

The data were corrected using the MC samples detailed in Section 4 for acceptance and the 
effects of detector response using the bin-by-bin method, in which the correction factor, as 
a function of an observable C in a given bin i, is Ci{0) = A^]''"^(C) /Nf^^i^O). The variable 
jyiiad^^-j jg number of events in the simulation passing the kinematic requirements on 
the hadronic final state described in Section 2 and Nf'^^{0) is the number of reconstructed 
events passing the selection requirements as detailed in Section 6. 

The results of a detailed analysis [15,39] of the possible sources of systematic uncertainty 
are listed below. Typical values for the systematic uncertainty are quoted for the cross 
sections as a function of x"^^, 

• varying the measured jet energies by ±1% [3,4,38] in the simulation, in accordance 
with the uncertainty in the jet energy scale, gave an uncertainty of =f5%; 

• the central correction factors were determined using the Pythia MC. The Herwig 
MC sample was used to assess the model dependency of this correction and gave an 
uncertainty of +4%, but up to -|-12% at lowest x"^^; 

• changing the values of the various cuts to remove backgrounds from DIS, cosmic-ray 
and beam-gas events gave a combined uncertainty of less than ±1%; 

• varying the fraction of direct processes between 34% and 70% of the total MC sample 
in order to describe each of the kinematic distributions gave an uncertainty of about 
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• changing the proton and photon PDFs to CTEQ4L [27] and WHIT2 [40] respectively 
in the MC samples gave an uncertainty of about ±1.5% and ±2.5%. 

The uncertainty in the cross sections due to the jet energy-scale uncertainty is correlated 
between bins and is therefore displayed separately as a shaded band in Figs. 4-13. All 
other systematic uncertainties were added in quadrature when displayed in these figures. 
The choice of MC sample also exhibited some correlation between bins and is hence given 
separately in Tables 2-20. In addition, an overall normalization uncertainty of 2.2% from 
the luminosity determination is not included in either the figures or tables. 



8 Results 



8.1 Dijet differential cross sections 

Differential cross-sections da/dET, da/dE^r^^^, da/dfj and da/dx°^^ are given in Tables 2- 
9 and shown in Figs. 4-7 for x"^^ above and below 0.75. For x°'^'^> 0.75, da/dEx and 
da/dE}^'^ fall by over three orders of magnitude over the Et and E^r^^'^ ranges measured 
and the jets are produced up to f/ ~ 2. For a;°'^^< 0.75, the slopes of da /dEx and da/ dE^^^^ 
are steeper, with the jets produced further forward in f/. It is interesting to note that in 
both regions of x^^, the data probe high values of x in the proton. 

The NLO QCD predictions, corrected for hadronization and using the AFG04 and CJK 
photon PDFs, are compared to the data. For x'^^> 0.75, the NLO QCD predictions 
describe the data well, although some differences in shape are observed for da/dEx and 
da/dE}^^^. Although measurements at high x'^^ are less sensitive to the structure of the 
photon, it is interesting to note that the prediction using the CJK photon PDF describes 
the Et spectrum somewhat better. The shapes for the f} and x'^^ distributions are also 
better reproduced using the CJK photon PDF. 

At low x"^*^, the difference in shapes between data and NLO QCD for da/dEx and 
da/dE^^^^ is more marked, as has been seen previously [4]. For the prediction using 
AFG04, the data and NLO agree in the lowest bin whereas the prediction is significantly 
lower at higher Et and E}^^^ . In contrast, the prediction from CJK is too high in the 
first bin, which dominates the cross section, but agrees well at higher Et and E}^^^ . For 
the fj and x'^^ distributions, the shapes are again better described by NLO QCD using 
the CJK photon PDF, although the normalization is too high. Sensitivity to the photon 
PDFs is discussed further in Section 8.4. 
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8.2 Measurement of dfr/d I A</>jj 



The cross-section da / d\/S.(j)^^\ is presented for 0:°'^'' above and below 0.75 in Tables 10 
and 11 and Fig. 8. For x°'^^> 0.75, the cross-section data fall by about three orders of 
magnitude in the cross section, more steeply than for x'^^< 0.75. The predictions from 
NLO QCD and also both Herwig and Pythia MC programs (plotted separately since 
the implementation of parton showers differs between the two programs) are compared to 
the data. The MC predictions are area normalized to the data in the measured kinematic 
region. At high x'^^, NLO QCD agrees with the data at highest jA^jjl, but it has a 
somewhat steeper fall off. The prediction from the Pythia MC program is similar to 
that for NLO QCD, whereas the prediction from the Herwig program describes the data 
well. For low the distribution for NLO QCD is much too steep and is significantly 
below the data for all values of |A(/)jj| except the highest bin. The prediction from the 
Pythia program is less steep, but still gives a poor description. The prediction from the 
Herwig program is in remarkable agreement with the data. 

The results and conclusions shown are qualitatively similar to those already seen in dijet 
photoproduction in which at least one of the jets was tagged as originating from a charm 
quark [11]. The results here confirm that the parton-shower model in Herwig gives a 
good simulation of high-order processes and suggest that a matching of it to NLO QCD 
would give a good description of the data in both shape and normalization. Should such 
a calculation or other high-order prediction become available, the distributions presented 
here would be ideal tests of their validity as they present inclusive quantities and also 
have higher precision compared to the previous result [1 1] . 



8.3 Optimized cross sections 

The cross-sections da/dx'^^, optimized to be most sensitive to the uncertainty on the 
gluon PDF in the proton, are given in Tables 12-19 and shown in Figs. 9 and 10 for x^^ 
above and below 0.75, respectively. The measurements cover a range in of about 0.1 
to 0.5. At high x'^^, the data are very well described by NLO QCD predictions. At low 
x'^^, the description by NLO QCD is poorer, particularly when using the AFG04 photon 
PDF. Generally the predictions with CJK describe the data better with the exception 
of the "Low- a;°^' 3" cross section. Inclusion of these high-x"^*^ data in future fits would 
constrain the proton PDFs further, in particular that of the gluon. To include the cross 
sections for low a systematic treatment of the photon PDFs and their uncertainty is 
needed. 
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8.4 Sensitivity to the photon PDFs 

As discussed in Section 8.1, the measured cross sections show sensitivity to the choice of 
photon PDFs. This is to be expected due to the extension further forward in pseudorapid- 
ity compared to previous measurements. This was investigated further, with the results 
presented in Figs. 11-13, where predictions with all five available parameterizations of 
the photon PDFs are compared to the data. In Table 20 and Fig. 11 the cross-section 
da/dx'^^ is shown. At high x'^^, all predictions are similar, as expected since there is little 
sensitivity to the photon structure in this region. Towards low x'^^, the predictions differ 
by up to 70%. The prediction from CJK deviates most from the other predictions and 
also from the data. The other predictions, although also exhibiting differences between 
each other of up to 25%, give a qualitatively similar description of the data. 

In Figs. 12 and 13, the cross-sections da/dE^, da/dx'^^ and da/df] are presented for 
x°^''<0.75, as shown previously in Figs. 4, 7 and 6, respectively, but here with additional 
predictions using different photon PDFs. For da/dEx, the prediction using CJK is much 
higher than the data in the first bin, but then agrees with the data for all subsequent bins. 
All photon PDFs have a similar shape, and none can reproduce the shape of the measured 
distribution. Apart from CJK, all PDFs are too low in the region 22.5 < Ex < 37.5 GeV. 
For the cross-section da/dx'^^, no prediction gives a satisfactory description of the data. 
The prediction from CJK is generally above the data by 20-30%, but describes the shape 
of the cross section reasonably well. All other predictions give a poor description of the 
shape, with cross sections which fall too rapidly to high For da/df], the prediction 
from CJK again gives the best description of the shape of the data, although it is too 
high in normalization. 

In summary, the data show a large sensitivity to the parameterization of the photon 
PDFs. The gluon PDF from CJK, in particular, differs from the others and this may give 
a hint of how to improve the photon PDFs. The data presented here should significantly 
improve the measurement of the gluon PDF of the photon, which is currently insufficiently 
constrained by the F2 data. 

9 Conclusions 

Dijet cross sections in photoproduction have been measured at high i?;^* and probe a wide 
range of x"^*" and The kinematic region is Q"^ < IGeV^, 142 < W^p < 293 GeV, 

^j.^" > 20 GeV, ^j?*^ > 15 GeV and -1 < r/j'^"'^ < 3, with at least one jet lying in the 
range between —1 and 2.5. In general, the data enriched in direct-photon events, at high 
x'^^, are well described by NLO QCD predictions. For the data enriched in resolved- 
photon events, at low x'^^, the data are less well described by NLO QCD predictions. 
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Predictions using different parameterizations of the photon parton density functions give 
a large spread in the region measured, with no parton density function giving an adequate 
description of the data. Therefore the data have the potential to improve the constraints 
on the parton densities in the proton and photon and should be used in future fits. The 
cross section in the difference of azimuthal angle of the two jets is intrinsically sensitive to 
high-order QCD processes and the data are poorly described by NLO QCD, particularly 
at low x*^^. Therefore the data should be compared with new calculations of higher orders, 
or simulations thereof. 
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Label 


x"*^^ cut 




cuts 


"High-x"^*^ 1" 
"Hieh-x°^'' 2" 

"High-x"^*^ 4" 
"Low-x°^" 1" 
"Low-x°^'^ 2" 
"Low-x°^'^ 3" 
"Low-x°^^ 4" 


x"^*^ > 0.75 
obs > 75 

x°^'^ > 0.75 
x"^*^ > 0.75 
x"^*^ < 0.75 
x°^'^ < 0.75 
x°^'^ < 0.75 
x°^'^ < 0.75 


< r/j^*^ < 1, 2 < r/j^*2 <- 3 

< r)j''*^ < 1 2 < 77^^=*^ < 3 

1 < ^^="'2 < 2 

-1 < r/j^" < 0, < < 1 
2 < r/j'^*! < 2.5, 2 < r/j^*2 < 3 
1 < < 2 

1 < < 2, 2 < r/j^*2 <- 3 
1 < < 2, 2 < <- 3 


^j.""'^ > 25, 15 GeV 
Ei?*^'^ > 20 15 GeV 
Ef^''^ > 30, 15 GeV 
^j?*^'^ > 20, 15 GeV 
^j?*^'^ > 20, 15 GeV 
Ei^'^'^ > 25, 15 GeV 
E}^'^'^ > 20, 15 GeV 
^j?"'^ > 25, 15 GeV 



Table 1: Kinematic regions of the "optimized" cross sections. 
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Err bin fGeV) 


daldE-r 


"Stat 




i^syst 








95 73 


+ 36 


+0.66 
-0.00 


+0.41 
-0.43 


+1.03 
-1.20 


n qc,c; -u n ni 7 
yj.cfoo —1— u.uxi 


22 5 27 5 


14.66 


± 0.28 


+0.00 
-0.28 


+0.42 
-0.26 


+0.60 
-0.65 


931 ± 008 


27.5, 32.5 


5.57 


± 0.18 


+0.09 
-0.00 


+0.14 
-0.24 


+0.30 
-0.19 


0.937 ± 0.029 


32.5, 37.5 


2.37 


± 0.12 


+0.00 
-0.03 


+0.15 
-0.04 


+0.11 
-0.11 


0.927 ± 0.012 


37.5, 42.5 


0.96 


± 0.07 


+0.02 
-0.00 


+0.06 
-0.03 


+0.07 
-0.03 


0.907 ± 0.034 


42.5, 55.5 


0.300 


± 0.024 


+0.000 
-0.004 


+0.004 
-0.018 


+0.016 
-0.020 


0.932 ± 0.044 


55.5, 70.5 


0.046 


± 0.009 


+0.006 
-0.000 


+0.001 
-0.003 


+0.003 
-0.003 


0.926 ± 0.029 


70.5, 90.5 


0.009 


± 0.004 


+0.001 
-0.000 


+0.001 
-0.002 


+0.000 
-0.002 


0.917 ± 0.085 



Table 2: Measured cross-section da/dEx for x'^^> 0.75. The statistical, (Jstat? MC 
model, 6mc, uncorrelated systematic, (5syst; and jet energy scale, 6es, uncertainties 
are shown separately. The hadronization correction factor, Chad? applied to the 
NLO QCD prediction is shown in the last column, where its uncertainty is half the 
spread between the values obtained using the Herwig and Pythia models. 



Et bin (GeV) 


da/dE-r 


"^stat 


^MC 


"^syst 


Ses (pb/GeV) 


Chad 


17.5, 22.5 


27.10 


± 0.36 


+0.49 
-0.00 


+0.18 
-1.31 


+1.45 

-1.42 


1.082 ± 0.045 


22.5, 27.5 


11.97 


± 0.24 


+0.07 
-0.00 


+0.21 
-0.66 


+0.56 
-0.74 


1.047 ± 0.009 


27.5, 32.5 


3.69 


± 0.14 


+0.17 
-0.00 


+0.10 
-0.23 


+0.27 
-0.18 


1.057 ± 0.016 


32.5, 37.5 


1.24 


± 0.08 


+0.03 
-0.00 


+0.06 
-0.12 


+0.07 
-0.09 


1.004 ± 0.024 


37.5, 42.5 


0.46 


± 0.05 


+0.03 
-0.00 


+0.01 
-0.05 


+0.04 
-0.03 


1.069 ± 0.043 


42.5, 55.5 


0.090 


± 0.013 


+0.005 
-0.000 


+0.009 
-0.010 


+0.008 
-0.007 


1.019 ± 0.015 


55.5, 70.5 


0.011 


± 0.005 


+0.004 
-0.000 


+0.006 
-0.002 


+0.001 
-0.001 


0.974 ± 0.064 



Table 3: Measured cross-section da/dEx for x'^^< 0.75. For further details, see 
the caption to Table 2. 
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^j?" bin (GeV) 


da 1 dE}^ 


^stat 








(pb/GeV) 




20, 26 


27.24 


± 0.33 


+0.18 
-0.00 


+0.56 
-0.54 


+1.05 
-1.22 




0.957 ± 0.021 


26, 32 


9.21 


± 0.20 


+0.17 
-0.00 


+0.21 
-0.15 


+0.49 
-0.37 




0.920 ± 0.011 


32, 38 


3.34 


± 0.12 


+0.00 
-0.05 


+0.16 
-0.12 


+0.14 
-0.17 




0.916 ± 0.024 


38, 44 


1.25 


± 0.07 


+0.03 
-0.00 


+0.15 
-0.03 


+0.07 
-0.06 




0.943 ± 0.005 


44, 55 


0.37 


± 0.03 


+0.00 
-0.00 


+0.01 
-0.03 


+0.02 
-0.03 




0.921 ± 0.035 


55, 70 


0.056 


± 0.009 


+0.008 
-0.000 


+0.004 
-0.003 


+0.007 
-0.002 




0.889 ± 0.051 


70, 90 


0.010 


± 0.004 


+0.004 
-0.000 


+0.004 
-0.001 


+0.002 
-0.000 




0.85 ± 0.11 



Table 4: Measured cross-section da / dE-)^ for x'J^^> 0.75. For further details, 
see the caption to Table 2. 



^j?" bin (GeV) 



da/dE^r^^^ 6s 



^MC 


^syst 


^ES 


+0.43 


+0.21 


+1.32 


-0.00 


-1.33 


-1.34 


+0.21 


+0.10 


+0.49 


-0.00 


-0.41 


-0.47 


+0.06 


+0.10 


+0.14 


-0.00 


-0.17 


-0.15 


+0.00 


+0.02 


+0.04 


-0.01 


-0.05 


-0.05 


+0.02 


+0.01 


+0.02 


-0.00 


-0.02 


-0.01 


+0.001 


+0.004 


+0.001 


-0.000 


-0.003 


-0.002 



(pb/GeV) 



had 



20, 26 
26, 32 
32, 38 
38, 44 
44, 55 
55, 70 



25.59 
8.11 
2.39 
0.72 
0.18 
0.018 



± 0.31 
± 0.18 
± 0.10 
± 0.05 
± 0.02 
± 0.006 



1.081 ± 0.043 
1.041 ± 0.015 
1.017 ± 0.025 
0.997 ± 0.006 
0.963 ± 0.027 
0.927 ± 0.033 



Table 5: Measured cross-section da/ dE^r^ for x^^< 0.75. For further details, 
see the caption to Table 2. 



7] bin 


da/drj 


^stat 


^MC 


^syst 


^ES (pb) 


Chad 


-0.50, 0.00 


4.8 


± 1.2 


+0.2 
-0.0 


+0.7 

-1.4 


+0.7 

-1.6 


0.551 ± 0.037 


0.00, 0.50 


90.1 


± 2.3 


+5.1 
-0.0 


+4.0 
-1.2 


+6.8 
-5.3 


0.892 ± 0.018 


0.50, 1.00 


177.8 


± 2.9 


+2.5 
-0.0 


+2.6 
-3.6 


+7.1 
-8.9 


0.940 ± 0.001 


1.00, 1.50 


167.6 


± 2.6 


+0.0 
-1.2 


+6.5 
-3.1 


+6.6 
-6.5 


0.952 ± 0.014 


1.50, 2.00 


59.0 


± 1.5 


+0.6 
-0.0 


+0.7 
-0.6 


+1.4 
-1.5 


1.079 ± 0.035 


2.00, 2.50 


2.8 


± 0.5 


+0.0 
-0.2 


+0.1 
-0.3 


+0.0 
-0.0 


1.062 ± 0.064 



Table 6: Measured cross-section da/drj for x?^^> 0.75. For further details, see 
the caption to Table 2. 
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7] bin 


da/drj 


<^stat 






^ES (pb) 


Chad 


0.00 


0.50 


7.2 


± 0.8 


+0.0 
-0.1 


+0.7 
-0.9 


+0.9 
-0.8 


1.052 ± 0.080 


0.50 


1.00 


65.9 


± 1.9 


+0.0 
-0.0 


+1.5 
-5.1 


+4.1 
-5.1 


1.074 ± 0.054 


1.00 


1.50 


144.0 


± 2.6 


+3.2 
-0.0 


+1.7 
-7.6 


+7.6 

-8.1 


1.080 ± 0.021 


1.50 


2.00 


146.8 


± 2.4 


+1.6 
-0.0 


+2.2 
-7.8 


+7.2 
-7.2 


1.063 ± 0.019 


2.00 


2.50 


71.3 


± 1.7 


+5.1 
-0.0 


+2.2 
-2.5 


+4.0 
-2.9 


1.062 ± 0.022 


2.50 


2.75 


18.4 


± 1.5 


+0.7 
-0.0 


+0.3 
-2.6 


+0.4 
-1.5 


1.066 ± 0.002 



Table 7: Measured cross-section da/ dt] for x'^'^< 0.75. For further details, see 
the caption to Table 2. 



bin 


da/dxf'' 


"^stat 




•^syst 


^ES (pb) 


Chad 


0.00, 0.05 


1260 


± 26 


+57 
-0 


+21 
-23 


+69 
-72 





902 


± 





025 


0.05, 0.10 


1960 


± 30 


+7 
-0 


+35 
-48 


+81 
-82 





932 


± 





007 


0.10, 0.15 


925 


± 20 


+0 

-1 


+60 
-12 


+27 
-41 





996 


± 





024 


0.15, 0.20 


468 


± 15 


+0 
-9 


+13 
-7 


+24 
-17 





999 


± 





015 


0.20, 0.25 


220 


± 11 


+0 
-4 


+12 
-5 


+6 
-9 





982 


± 





012 


0.25, 0.30 


104.9 


± 8.4 


+0.0 
-1.3 


+2.9 
-10.8 


+5.1 
-4.1 





963 


± 





015 


0.30, 0.35 


45.0 


± 5.6 


+1.5 
-0.0 


+3.4 
-1.0 


+2.4 
-1.2 


1 


063 


± 





023 


0.35, 0.40 


23.2 


± 4.1 


+0.0 
-0.9 


+0.5 
-0.9 


+0.6 

-1.6 


1 


027 


± 





008 


0.40, 0.45 


8.7 


± 2.4 


+0.9 
-0.0 


+4.0 
-0.5 


+1.0 
-0.1 


1 


010 


± 





020 


0.45, 0.50 


3.2 


± 1.4 


+0.0 
-0.3 


+2.5 
-1.0 


+0.2 
-0.2 


1 


006 


± 





016 


0.50, 1.00 


0.40 


± 0.17 


+0.08 
-0.00 


+0.08 
-0.21 


+0.06 
-0.01 





987 


± 





018 



Table 8: Measured cross-section da/ dx'^^ for x'^^> 0.75. For further details, see 
the caption to Table 2. 
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™obs 

0.00 
0.05 
0.10 
0.15 
0.20 
0.25 
0.30 
0.35 
0.40 
0.45 
0.50 



bin 

0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

1.00 



da/dx°J 
236 
1131 
1120 
829 
581 
302 
146.8 
65.5 
24.6 
9.6 
0.86 



obs 



^stat 
± 12 

± 24 
± 22 
± 19 
± 17 
± 12 
± 9.4 
± 6.6 
± 4.2 
± 2.7 
± 0.27 





^syst 


^ES (pb) 


Chad 


+2 
-0 


+17 
-24 


+18 
-19 


1 


103 


± 





092 


+0 
-0 


+19 
-76 


+55 
-70 


1 


063 


± 





046 


+19 

-0 


+37 
-63 


+56 
-61 


1 


086 


± 





022 


+12 
-0 


+7 
-37 


+46 
-37 


1 


074 


± 





001 


+14 

-0 


+5 
-49 


+31 
-30 


1 


053 


± 





001 


+31 
-0 


+25 
-10 


+17 
-13 


1 


052 


± 





052 


+8.3 
-0.0 


+4.2 
-6.2 


+7.0 
-9.7 


1 


052 


± 





014 


+0.0 

-0.3 


+0.6 
-15.0 


+3.9 
-4.2 


1 


041 


± 





008 


+1.1 
-0.0 


+4.8 
-2.2 


+0.4 
-3.0 


1 


036 


± 





004 


+0.0 

-0.7 


+0.7 
-2.3 


+1.7 
-0.2 


1 


020 


± 





005 


+0.09 
-0.00 


+0.32 
-0.09 


+0.07 
-0.10 


1 


012 


± 





034 



Table 9: Measured cross 
the caption to Table 2. 



■section da/dx°^^ for x'^^< 0. 



75. For further details, see 



|A</>ij| bin 


da/d\A(l^^\ 


^stat 




^syst 


^ES 


(pb/rad) 


Chad 


1.83, 2.09 


1.7 


± 0.5 


+0.1 
-0.0 


+0.2 
-0.5 


+0.1 
-0.2 




0.65 ± 0.11 


2.09, 2.36 


7.8 


± 1.0 


+0.0 
-0.0 


+1.2 
-0.6 


+0.6 
-0.6 




0.729 ± 0.059 


2.36, 2.62 


36.1 


± 2.2 


+0.2 
-0.0 


+1.6 
-1.7 


+2.1 
-1.8 




0.826 ± 0.013 


2.62, 2.88 


132.9 


± 3.9 


+5.8 
-0.0 


+5.9 
-2.7 


+6.6 
-8.3 




0.868 ± 0.008 


2.88, 3.14 


779.1 


± 8.1 


+4.0 
-0.0 


+15.0 
-13.3 


+31.8 
-33.6 




0.984 ± 0.015 



Table 10: Measured cross-section da/ d\A(l)^^\ for x'^^> 0.75. For further details, 
see the caption to Table 2. 



|A0jj 


bin 


d(r/d\A^^\ 


^stat 


0.00, 


1.57 


0.26 


± 0.07 


1.57, 


1.83 


2.9 


± 0.6 


1.83, 


2.09 


6.6 


± 0.8 


2.09, 


2.36 


28.2 


± 1.7 


2.36, 


2.62 


78.4 


± 2.8 


2.62, 


2.88 


203.2 


± 4.5 


2.88, 


3.14 


528.6 


± 6.7 





^syst 




(pb/rad) 


Chad 


+0.05 
-0.00 


+0.02 
-0.02 


+0.04 
-0.02 




0.84 ± 





15 


+0.3 
-0.0 


+0.6 
-0.1 


+0.1 
-0.3 




0.869 ± 





083 


+0.2 
-0.0 


+0.4 
-0.2 


+0.3 
-0.6 




0.910 ± 





031 


+0.0 
-0.5 


+0.6 
-2.3 


+2.4 
-1.3 




0.959 ± 





004 


+1.2 
-0.0 


+3.5 
-1.0 


+4.3 
-5.3 




0.988 ± 





006 


+0.0 
-1.1 


+0.6 
-8.6 


+10.4 
-13.4 




1.006 ± 





015 


+16.5 
-0.0 


+6.0 
-36.5 


+28.1 
-26.4 




1.069 ± 





020 



Table 11: Measured cross-section da/ d\A(f^^\ for x'^^< 0.75. For further details, 
see the caption to Table 2. 
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™obs 


bin 


da/dxf^ 


f^stat 




^syst 




(pb) 


Chad 


0.1, 


0.2 


80.9 


± 4.2 


+0.0 

-3.4 


+3.8 
-6.1 


+3.8 
-3.4 




0.957 ± 0.010 


0.2, 


0.3 


51.6 


± 3.5 


+0.0 
-1.0 


+3.1 
-2.0 


+2.4 
-2.1 




0.974 ± 0.059 


0.3, 


0.4 


12.6 


± 2.1 


+0.0 
-0.0 


+1.0 
-0.9 


+0.6 
-0.9 




0.962 ± 0.010 


0.4, 


0.5 


2.1 


± 1.0 


+1.0 
-0.0 


+1.0 
-0.3 


+0.2 
-0.1 




0.953 ± 0.024 



Table 12: Measured cross-section da/dzf"" for x°^''> 0.75 ("High-xf'' 1"). For 
further details, see the caption to Table 2. 



xf''' bin 


da/dxf^ 


^stat 




^syst 


^ES 


(pb) 


Chad 


0.0, 0.1 


10.1 


± 1.6 


+0.1 
-0.0 


+0.6 

-0.5 


+0.7 
-0.2 




0.961 ± 0.037 


0.1, 0.2 


238.9 


± 7.1 


+0.0 
-5.2 


+15.0 
-6.8 


+9.7 
-10.8 




1.006 ± 0.021 


0.2, 0.3 


77.0 


± 4.5 


+0.0 
-2.4 


+6.7 
-1.9 


+3.6 
-2.7 




1.005 ± 0.026 


0.3, 0.4 


12.6 


± 2.1 


+0.0 
-0.0 


+0.9 
-0.9 


+0.6 
-0.9 




0.964 ± 0.009 


0.4, 0.5 


2.1 


± 1.0 


+ 1.0 
-0.0 


+1.0 
-0.3 


+0.2 
-0.1 




0.953 ± 0.024 



Table 13: Measured cross-section da/dxf"" for x°^''> 0.75 ("High-x"^'' 2"). For 
further details, see the caption to Table 2. 



x°p^ bin 


da/dxf^ 


f^stat 




^syst 




(pb) 


Chad 


0.0, 0.1 


2.1 


± 0.8 


+0.4 
-0.0 


+1.4 
-0.1 


+0.1 
-0.1 




0.914 ± 0.014 


0.1, 0.2 


55.9 


± 3.5 


+0.1 

-0.0 


+1.2 
-2.7 


+2.3 
-1.4 




0.974 ± 0.006 


0.2, 0.3 


20.5 


± 2.1 


+0.9 
-0.0 


+0.3 
-3.0 


+0.7 
-0.8 




0.988 ± 0.011 


0.3, 0.4 


2.4 


± 0.7 


+0.0 
-0.0 


+0.1 
-0.4 


+0.1 
-0.1 




1.007 ± 0.046 



Table 14: Measured cross-section da/dxf"" for x°^''> 0.75 ("High-x"^'' 3"). For 
further details, see the caption to Table 2. 



x°^^ bin 


da/dxf"" 4tat ^Mc 4yst 5es (pb) 


Chad 


0.0, 0.1 


1Q8n -1-8 8 +10-9 +2-9 +18.7 
iy».U ± «.« „o.O -2.3 -16.0 


0.832 ± 0.017 



Table 15: Measured cross-section da/dx^^^ for x°^^> 0.75 ("High-x°^^ 4")- Por 
further details, see the caption to Table 2. 
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Xp^^ bin 


da/clxf"' 


^stat 




'^syst 




(pb) 


Chad 


0.1, 0.2 


15.0 


± 2.0 


+0.8 
-0.0 


+2.2 
-0.5 


+0.5 
-0.3 




1.004 ± 0.099 


0.2, 0.3 


89.4 


± 4.6 


+13.4 
-0.0 


+1.5 
-4.1 


+4.3 
-3.9 




1.030 ± 0.003 


0.3, 0.4 


46.7 


± 3.8 


+2.3 
-0.0 


+0.4 
-4.3 


+1.8 
-3.3 




1.070 ± 0.090 


0.4, 0.5 


7.0 


± 1.5 


+0.4 
-0.0 


+0.2 
-0.6 


+0.1 
-0.9 




0.960 ± 0.083 


0.5, 1.0 


0.48 


± 0.20 


+0.00 
-0.04 


+0.04 
-0.09 


+0.03 
-0.05 




1.024 ± 0.027 



Table 16: Measured cross-section da/dx^^ for x^^< 0.75 ("Low-x^^ 1"). For 
further details, see the caption to Table 2. 



xf"" bin 


da/dxf^ 


^stat 






^ES (pb) 


Chad 


0.0, 0.1 


19.5 


± 2.3 


+1.5 
-0.0 


+0.8 
-3.0 


+0.4 
-1.8 


0.876 ± 0.076 


0.1, 0.2 


117.6 


± 5.0 


+2.0 
-0.0 


+4.7 
-9.7 


+5.5 
-5.3 


1.048 ± 0.014 


0.2, 0.3 


12.6 


± 1.7 


+0.6 
-0.0 


+0.6 
-1.9 


+0.7 
-0.7 


1.116 ± 0.085 



Table 17: Measured cross-section da/dxf'^ for x'^^< 0.75 ("Low-x^^ 2"). For 
further details, see the caption to Table 2. 



x°p^ bin 


da/dxf^ 


^stat 




'^syst 


(^ES (pb) 


Chad 


0.1, 0.2 


278.4 


± 7.6 


+4.2 
-0.0 


+4.6 
-12.7 


+13.5 
-12.4 


1.087 ± 0.015 


0.2, 0.3 


235.2 


± 7.1 


+10.3 
-0.0 


+2.1 
-9.6 


+12.2 
-10.3 


1.077 ± 0.030 


0.3, 0.4 


47.8 


± 3.6 


+0.7 
-0.0 


+0.8 
-3.4 


+2.8 
-2.6 


0.999 ± 0.064 


0.4, 0.5 


8.3 


± 1.6 


+0.0 
-0.1 


+1.7 
-0.6 


+0.7 
-0.6 


1.037 ± 0.020 


0.5, 1.0 


0.28 


± 0.14 


+0.15 
-0.0 


+0.19 
-0.04 


+0.07 
-0.01 


1.003 ± 0.037 



Table 18: Measured cross-section da/dxf^ for xf''< 0.75 ("Low-xf^ 3"). For 
further details, see the caption to Table 2. 



x'^^ bin 


da/dxf^ 


f^stat 


0.1, 0.2 


71.3 


± 4.1 


0.2, 0.3 


120.4 


± 5.0 


0.3, 0.4 


45.0 


± 3.4 


0.4, 0.5 


8.3 


± 1.6 


0.5, 1.0 


0.28 


± 0.14 



+1.8 
-0.0 
+5.6 
-0.0 
+0.3 
-0.0 
+0.0 
-0.1 
+0.15 
-0.00 



"syst 

+2.6 
-4.6 
+2.6 
-6.3 
+1.9 
-3.3 
+1.7 
-0.6 
+0.19 
-0.04 



Pes 

+4.2 
-3.4 

+7.3 
-4.6 
+1.8 
-3.2 
+0.7 
-0.6 
+0.07 
-0.01 



(pb) 



Chad 

1.066 ± 0.052 
1.042 ± 0.021 
1.013 ± 0.059 
1.037 ± 0.020 
1.003 ± 0.037 



Table 19: Measured cross-section da/dx°^^ for x°^^< 0.75 ("Low-x'^^ 4")- For 
further details, see the caption to Table 2. 
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x'^^^ bin 


da/dx"^^' 


^stat 






^ES (pb) 


Chad 


0.1, 0.2 


169.5 


± 6.8 


+19.6 
-0.0 


+2.3 
-7.4 


+14.7 
-12.6 


1.081 ± 0.046 


0.2, 0.3 


271.6 


± 8.0 


+12.0 
-0.0 


+1.7 
-8.2 


+17.1 
-14.3 


1.042 ± 0.056 


0.3, 0.4 


325.7 


± 8.9 


+0.3 
-0.0 


+2.5 

-15.2 


+16.2 
-16.3 


1.065 ± 0.017 


0.4, 0.5 


346.6 


± 9.3 


+7.2 
-0.0 


+7.6 
-15.3 


+17.2 
-19.0 


1.058 ± 0.023 


0.5, 0.6 


385 


± 10 


+3 
-0 


+4 
-21 


+20 
-19 


1.072 ± 0.016 


0.6, 0.7 


458 


± 11 


+3 

-0 


+17 
-30 


+20 

-24 


1.089 ± 0.028 


0.7, 0.8 


557 


± 12 


+1 
-0 


+16 
-55 


+28 
-29 


1.087 ± 0.011 


0.8, 1.0 


1106 


± 11 


+15 

-0 


+32 
-21 


+47 
-48 


0.940 ± 0.018 



Table 20: Measured cross-section da/dx'^^. For further details, see the caption 
to Table 2. 
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Figure 1: Examples of (a) direct and (b) resolved dijet photoproduction diagrams 
in positron-proton, ep, collisions in LO QCD. This direct-photon process is the 
collision of a photon, 7, and gluon, g from the proton. This resolved-photon process 
is a collision of a parton from the photon and a gluon, g, from the proton. 
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Figure 2: The theoretical uncertainties (see Section 5) for sample distributions: 
(a) (h) Et forx°^^< 0.75, (c) "Low-x°^' 3" and (d) "High-x"^' 2", which are 
defined in Table 1. The uncertainties are the total (outer shaded hand), that from 
varying fi^ (inner shaded band), the experimental uncertainties of data propagated 
in the ZEUS-JETS fit (solid lines) and using the most different photon PDF, CJK 
(dashed line) instead of AFG04- 
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Figure 3: Predictions of the fraction of the cross section initiated by gluons for 



sample distributions: (a) x'^^ , (b) Et for x"^^ 



< 0.75, (c) "Low-xf'' 3" and (d) 



"High-x'^^ 2", which are defined in Table 1. The gluon fractions are from the 
proton using the CTEQ5M1 PDF (long-dashed line), and from the photon using 
the AFGO4 (solid line) and CJK PDFs (short-dashed line). 
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Figure 4: Measured cross-section da/dEx for (a) x'^^> 0.75 and (b) x'^^< 
0.75 compared with NLO QCD predictions using the AFGO4 (solid line) and CJK 
(dashed line) photon PDFs. The data (dots) are shown with statistical (inner bars) 
and statistical and systematic uncertainties added in quadrature ( outer bars ) along 
with the jet energy- scale (Jet ES) uncertainty (shaded band). The NLO QCD pre- 
dictions are shown (NLO QCD ® HAD) multiplied by the hadronization correc- 
tions, Chad; discussed in Section 5. The predictions using AFCO4 are also shown 
with their associated uncertainties (shaded histogram) as discussed in Section 5. 
The ratios to the prediction using the AFCO4 photon PDF are shown at the bottom 
of the figure. 
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Figure 5: Measured cross-section da/dE}^''^ for (a) x°^^> 0.75 and (h) 
0.75. For further details, see the caption to Fig. 4- 
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Figure 6: Measured cross-section da/drj for (a) x'^^> 0.75 and (h) x'^^< 0.75. 
For further details, see the caption to Fig. 4- 
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Figure 7: Measured cross-section da/dx^^^ for (a) x'^^> 0.75 and (h) x°^^< 0.75. 
For further details, see the caption to Fig. 4- 
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Figure 8: Measured cross-section da/ dA\(f)^^\ for (a) x'^^> 0.75 and (h) x^^< 0.75 
compared with NLO QCD predictions using the AFG04 (solid line) photon PDF. 
Predictions from the MC programs Herwig (dot-dashed) and Pythia (dashed), 
area normalized to the data by the factors given, are also shown. The data (dots) 
are shown with statistical (inner bars) and statistical and systematic uncertainties 
added in quadrature (outer bars) along with the jet energy-scale (Jet ES) uncer- 
tainty (shaded band). The NLO QCD predictions are shown (NLO QCD ® HAD) 
multiplied by the hadronization corrections, Chad? discussed in Section 5. The pre- 
dictions using AFGO4 are also shown with their associated uncertainties (shaded 
histogram) as discussed in Section 5. 
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Figure 9: Optimized cross-sections da/dx'^^ for x'^^> 0.75 in the kinematic 
regions defined in Table 1. For further details, see the caption to Fig. 4- 
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Figure 10: Optimized cross-sections da/dx'^^ for x'^^< 0.75 in the kinematic 
regions defined in Table 1. For further details, see the caption to Fig. 4- 



33 



ZEUS 



^1400 
n 

o ^- 

-giooo 

■D 800 



1 — I — \ — \ — I — I — \ — \ — r 

• ZEUS 82 pb"^ 
NLO® HAD: 
— AFG04 
■-- CJK 

AFG 

h-' - GRV 
SAL 

■ Jet ES uncertainty 




400 
200 



rObs 

7 



Figure 11: Measured cross-section for da/dx" compared with NLO QCD pre- 
dictions using the AFGO4 (solid line), CJK (dashed line), AFG (dotted line), GRV 
(dashed and double-dotted line) and SAL (dashed and single-dotted line) photon 
PDFs. The data (dots) are shown with statistical (inner bars) and statistical and 
systematic uncertainties added in quadrature (outer bars) along with the jet energy- 
scale (Jet ES) uncertainty (shaded band). The NLO QCD predictions are shown 
(NLO QCD ® HAD) multiplied by the hadronization corrections, Chad? discussed in 
Section 5. The predictions using AFGO4 are also shown with their associated uncer- 
tainties (shaded histogram) as discussed in Section 5. The ratios to the prediction 

using the AFGO4 photon PDF are shown at the bottom of the figure. 
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Figure 12: Measured cross-section for da/ (IEt for x'^^< 0.75. For further details, 
see the caption to Fig. 11. 
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Figure 13: Measured cross-section for (a) da/dx^^ and (b) da/drj both for x'^^< 
0.75. For further details, see the caption to Fig. 11. 
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